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q-i General relativity allows a variety of future singularities to occur 

in the evolution of the universe. At these future singularities, the 
universe will end in a singular state after a finite proper time and 
geometrical invariants of the space time will diverge. One question 
. that naturally arises with respect to these cosmological scenarios is 

1 the following: can quantum effects lead to the avoidance of these 

future singularities? We analyze this problem considering massless 
C*"} ■ and conformally coupled scalar fields in an isotropic and homoge- 

neous background leading to future singularities. It is shown that 
near strong, big rip-type singularities, with violation of the energy 
conditions, the quantum effects are very important, while near some 
milder classes of singularity like the sudden singularity, which pre- 
serve the energy conditions, quantum effects are irrelevant. 
PACS number: 98.80-k 
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1 Introduction 

It is general believed that today the Universe is in a stage of accelerated ex- 
pansion. The primary evidence for this accelerated phase of cosmic evolution 
came from the use of the supernova type la as standard candles to measure 
distances in the universe. Supernova are very bright objects that can be seen 
to great distances. Presently, measurements of supernova type la up to z ~ 1.8 
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are available. The fact that these very distant supernova appear dimmer than 
would be expected in a pure matter-dominated universe led to the conclusion 
that the universe is accelerating. Much speculation has arisen as to the source of 
this acceleration and even its reality has been questioned. But, the fact that the 
spectrum of anisotropy of the cosmic microwave background radiation indicates 
that the spatial section of the geometry of the universe must be almost flat, 
while the observations of the virialized system (galaxies, clusters of galaxies, 
etc.) indicate a low density of the universe, implies indirectly the acceleration 
of the universe. In fact, in order to complete the cosmic energy budget and 
have flat spatial sections, it is necessary to include a component that does not 
agglomerate locally, but remains as a smooth component of universe. To have 
this feature, this component must have negative pressure, and consequently it 
must dominate the matter content of the universe asymptotically, driving the 
accelerated expansion in the later phases of the cosmic evolution. 

In brief, to explain the acceleration of the universe, an exotic component in 
the cosmic budget exhibiting negative pressure is needed. This exotic compo- 
nent is named dark energy. We must remember that observations require also 
a second non-baryonic component, called dark matter, with zero pressure, nec- 
essary to explain conveniently the formation of structures in the universe and 
the dynamics of local, virialized system, like galaxies and clusters of galaxies. 
However, while there exists a lot of reasonable candidates to represent dark 
matter (neutralinos, axions, sterile neutrinos, primordial black holes etc - for 
a review, see reference [I]), it is not clear what kind of fluid or field would 
constitute dark energy. The first natural candidate to be evoked has been the 
cosmological constant, seen as a phenomenological manifestation of the vacuum 
energy of quantum fields existing in the universe. However, this very attrac- 
tive possibility is plagued with fine -tunings and poorly understood issues, (see 
the classical references [2]). In spite of this, the cosmological constant remains 
the most popular candidate to represent dark energy, leading to the so-called 
ACDM model, highly supported by observations. A simple extension of general 
relativity, created by restricting the variational principle for deriving the Ein- 
stein equations to causal variations, leads to a general prediction that there is a 
cosmological constant with the observed value and a prediction that the value 
of the dimensionless curvature will be observed to be -0.0055 [3]. 

The cosmological constant implies an equation of state such that a x — 
Px/px = — 1, where p x is the density and p x is the pressure (the subscript x 
denotes a dark energy component). Observations will never give an exactly 
precise result for the equation of state of dark energy (or any component of the 
universe), due to statistical and systematic errors. But, if these dispersions in 
the evaluation of the dark energy equation of state can be reduced substantially, 
a clear case for the simple cosmological constant can emerge. We are far from 
this situation. But, in the search for the determination of a x some curious re- 
sults have appeared: while the estimations lead to values near a x = — 1, more 
negative values are highly admitted. This fact has raised a dramatic specula- 
tion: could dark energy have an equation of state such that a x < — 1? In fact, 
for example, if a x is constant and the spatial sections of the universe is flat, the 
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recent results of 7- years WMAP observations indicate that a x = — 1.10±0.14 
0] at la. 

If a x < — 1, all energy conditions are violated, in particular the null energy 
condition, which requires p x + p x > 0. From the conservation of the energy- 
momentum tensor in an expanding universe, we obtain the equation, 

px + Z-ipx+Px) = 0, (1) 

a 

a where a is the scale factor of the universe, and dots mean derivatives with 
respect to the cosmic time. From this expression, it comes out that the violation 
of the null energy condition implies that the energy density of the dark energy 
fluid grows with the expansion of the universe, instead of decreasing. A remark- 
able consequence of such behaviour is that any universe dominated by such an 
exotic fluid will inevitably hit a singularity in its future evolution, after a finite 
proper time. At this future singularity, the expansion scale factor a and the en- 
ergy density p x both diverge. This highly singular state in the future evolution 
of the universe has been called big rip singularity [5] , and it is a remarkable and 
plausible example of what are called future cosmological singularities. 

However, the situation is more complicated, and future singularities can exist 
even if the energy conditions are not violated. The possibility of future singu- 
larities was pointed out for the first time in reference [5] and their occurrence 
does not necessarily require the violation of the energy conditions. There is also 
a later discussion in reference [7] . A recent example of future singularities which 
does not violate the null energy condition is the 'big brake' singularity, which 
emerges from the DBI action [8] [9] . The big brake singularity has the curious 
property that it can be traversed by a pointlike particle [lOj . A milder type 
of finite time singularity is the so-called " sudden singularity" [TT] . The sudden 
future singularity occurs without violating any energy condition (so p + p > 
and p + 3p > at all times) . This singularity is characterized by a finite value 
for the scale factor, its first time derivative, and the density, while the second 
derivative of the scale factor and the pressure diverge at finite time [TT]. They 
are singularities of the weak sort discussed by Tipler [T2] and Krolak [T3] . They 
cannot occur if dp/ dp is continuous. 

In general, it is believed that the fate of the universe near any singularity 
(past or future) must be affected by quantum effects arising in the extreme 
conditions that exist in its spacetime neighborhood. In the case of the big rip, 
this problem has been treated for example, in references [HI [XSJ [16]. In these 
investigations, it was found that the quantum effects are important. But, the 
conclusions concerning the back reaction of the quantum effects on the evolution 
of the universe were harder to decide unambiguously. In the case of the sudden 
singularity, quantum effects were studied in references [TT] [T5] [Jj5], and the 
results indicated that quantum effects do not change the evolution near the 
singularity. These results were obtained for massive and massless scalar fields. 

Here, the fate of quantum effects near future singularities is reviewed. We 
begin, in next section, by describing the some general class of future singularities. 
In section 3 the evolution of a scalar field in the background of big rip and sudden 
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future singularities is determined. In section 4, the structure of the background 
solution is shown, and the solution for the scalar field equation is found. In 
section 5, quantum effects near future singularities are studied. In section 6 we 
present our conclusions. 

2 Future singularities 

Generally, a singularity in geometric theories of gravitation can be characterized 
by the divergence of some curvature invariants or, alternatively, by the incom- 
pleteness of the geodesic trajectories. Cosmology is commonly studied using a 
homogeneous and isotropic space-time. In this case, all dynamics is encoded 
in the behavior of a single function, the scale factor a(t), determined by some 
matter density p{t) endowed with a pressure p(t). In this case, singularities 
appear as a divergence in the Hubble factor H = - and in the density p(t). 
A traditional example is the initial big bang singularity, for which H(t) — > oo 
and p{t) — > oo as t — > 0. Of course, at the big bang singularity the curvature 
scalars diverge and the geodesies "begin" in the singularity, implying the the 
space-time is geodesically incomplete. 

The future singularities which occur at the end of a cosmological evolution, 
have in principle many of the features as the big bang singularity, including the 
divergence of the curvature invariants, which are displayed by many types of 
singularity. However, geodesic incompleteness is not a universal characteristic 
of these singularities and nor is a divergence in the Hubble parameter or in the 
density. The big rip singularity, which requires the violation of all energy con- 
ditions, bears a close resemblance to the big bang singularity, since, besides the 
divergence of curvature invariants, it is incomplete geodesically, and the energy 
density and the Hubble parameter both diverge. This is direct consequence of 
the violation of the null energy condition. As already remarked in the intro- 
duction, since p + p is negative when the null energy condition is violated, the 
energy density grows as the universe expands, in contrast to the usual situation. 
As consequence, R — > oo, H — > oo and p — > oo after a finite time t s . Moreover, 
the geodesies terminate at this singularity, and are inextendible. 

The case of the big rip singularity can be described in a simple way by 
assuming a fluid (named from now on phantom fluid) with an equation of state 
of the type 

p = ap, a < —1. (2) 
Solving the conservation equation, the energy density scales as 

pota- 3(1+a K (3) 

Since, a < —1, p grows as a increases. 

This singularity occurs in a finite proper time. This can be seen more simply 
in the following way. In terms of the cosmic time, the scale factor behaves as 

floct 5 * 1 ^). (4) 
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Since 1 + a is negative, this solution can represent an expanding universe if 
— oo < t < 0_. It must be remarked that at t — > —oo there is no singularity, as 
we will see below. If we live in a moment to < 0, the time to elapse from t = to 
to t = 0_ is obviously finite. If there is any other form of "normal" matter, the 
phantom fluid will always dominate asymptotically, since the density of normal 
forms of matter decreases with the expansion. In reference [I] a simple model 
including pressureless matter and phantom fluid has been described. On the 
other hand, since R cx ™, p oc h and H cx j, all these quantities diverge as 
t — » 0_. It is in this sense that the big rip singularity can be considered as the 
reverse of the big bang singularity, with a Minkowskian asymptotic space-time 
as a — > and a singular state as a — > oo. 

There are other classes of mild future singularities. We understand a "mild 
singularity" to be a singularity that exhibits a divergence in the curvature in- 
variants (which is a requirement), but perhaps with no divergence in the density. 
If the density does not diverge then Einstein's equations imply that the Hubble 
function also does not diverge at the singularity. At same time, since density is 
connected with the Hubble function, it must reach a finite value at singularity. 
However, the Ricci scalar is given by 
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In order to have a divergence in R, we see that d must diverge. Through 
Einstein's equations, this implies that the pressure must diverge. This curious 
structure may allow that the geodesies to be continued through the singularity. 

One example of these class of "mild future singularities" is the sudden sin- 
gularity, which can be described by the following expression for the scale fac- 
tor pen [is on]: 

a(t)=(i)V-l) + l-(l-i)", (6) 

where t s is the time where the sudden singularity occurs, and a s is the value of 
the scale factor at this moment. Moreover, < q < 1 and 1 < n < 2 where q 
and n are free constants and no specific relation is assumed between the pressure 
p and the density p. 

The sudden singularity does not require any violation of the energy condi- 
tions. Pressure can remain positive. The only requirement is the divergence of 
pressure in a given finite time. One difficult with such scenario is to find a type 
of matter that could satisfy these requirements. While this seems to be difficult 
using ordinary types of fluids, non-standard self-interacting scalar fields may 
fulfill the requirements needed for a sudden singularity scenario. 

In fact, a self-interacting scalar field which presents non-standard coupling 
structures may lead to new and unexpected scenarios. One example is given by 
the DBI action, 

£ = V=gV(T) VI - T. tP T'P, (7) 

where T is a tachyonic scalar field, and V(T) is a potential term. This action 
can emerge from some specific configuration of string theory. The choice of 
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the potential V(T) may lead to very new kind of cosmological scenario. One 
example has been given in [7], where a new kind of future singularity has been 
exhibited: the big brake singularity. In this singularity, the Hubble parameter 
vanishes on the singularity, while the second derivative of the scale factor goes 
to — oo, from which the name big brake. Geodesies can also traverse the big 
brake singularity. 

From now on, we will concentrate our analysis on the big rip and sudden 
singularities. 



3 The master equation 

We will consider the quantum creation of particles near a given future cosmo- 
logical singularity. Two types of future singularity will be analyzed: the big rip 
singularity and the sudden singularity. Let us initially consider a general mas- 
sive, non-minimally coupled scalar field tf>, giving by the following Lagrangian: 

C = - \m 2 4? + |i?0 2 . (8) 

The conformal coupling corresponds to £ = g . From the Lagrangian (jHJ we can 
deduce the following field equation: 

□0 + m 2 4> - £R0 = 0. (9) 

For a flat FLRW metric, this equation reduces to 



a I . I t.2 , 2 2 — a 



2 — f +^ r + mV + 6£— U = 0. (10) 



a 



In general, at future singularities the second derivative of the scale factor 
diverges, while the first one remains finite. Hence, it is convenient to transform 
equation (I10[) in order to eliminate this possibly singular term. This can be 
achieved by defining This transformation leads to the equation 

X " +2(l-6^- X ' 



fc 2 +mV + 6£(6f-l)f^ |x = 0. (11) 

The equation (fTTj) simplifies considerably for a conformal coupling, £ = g. If, 
a massless field is also considered, equation (fTTj) takes the form of an harmonic 
oscillator equation. 

Particle production near the big rip has been analyzed the references [T31 [TH 
115) . considering only massless scalar particle, while the corresponding analysis 
for the sudden singularity has been performed in reference |16) . In reference |17) 
the massive case with conformal coupling in the background of the sudden sin- 
gularity has been studied. We will review the results obtained in these references 
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later. We must remark that the non-minimal coupling should strengthen any 
quantum effects near the singularity, and also introduces some technical features 
that leads to exact solutions for the problem. 

Under variations with respect to the metric, the Lagrangian ([5]) gives the 
following momentum-energy tensor: 



where 



G^u — Rpv — -^9nvR (13) 

is the Einstein tensor. For a conformal coupling (£ = |), this expression reduces 
to, 

2 1 1 2 2 

r-g^n^+hs^tp. (14) 

In the case of minimally coupled massive field (£ = 0), the energy- momentum 
tensor reduces to 

TfMu = <j>;n<f>;v - \flpv$;p$ P + g M „ <f) 2 . (15) 

From now on, we will consider two main cases: m^O and £ = i for the 
massive conformally coupled scalar field, and m = and £ = 0, for the massless 
minimally coupled scalar field. 

4 The cosmological background and the solu- 
tions of the master equation 

The first task in analyzing the quantization of the scalar field equations is to 
determine the background where the quantum field evolves. As stated above, 
two cases will be considered: the big rip and the sudden singularity. 

4.1 The big rip background 

Let us consider the flat Friedmann equation and the conservation law for a fluid 
with an equation of state p = ap: 

a'\ 2 8ttG 2 

-a) = — pa > (16) 
p' + 3(l + a)p = 0. (17) 
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The general solution for the scale factor is given by a = ao\rj\ 1+3q (77 is the 
conformal time defined by dt — a(r})dr]), while the density behaves as p = 
poa~ 3 ( 1+a \ If a < — | then accelerated expansion occurs. Moreover, if a < 
— 1 (phantom regime) then the matter density grows during this accelerated 
expansion, leading to a big rip. One of the main aspects of these accelerating 
solutions is that the conformal time takes values in the interval — oo > 77 > 0_. 
Note that for 77 — > —00 the scale factor goes to zero, and the Minkowski space- 
time is asymptotically approached in the phantom regime, since p — > in this 
regime. This property will be important in order to fix the initial conditions. 

4.2 Sudden singularity: two cosmological eras 

For the sudden singularity, the situation is more complex due to the complicated 
form of the scale factor. But, the scale factor ^ admits two asymptotic forms, 
which we will call the primordial phase and the singular phase, respectively. 



Primordial phase, t — > 0: 



a -> ( -) (o.-l), (18) 



q ( t x 



(a s - 1), (19) 



9 , ,Jt 



• Singular phase, t — > t s 



q-2 

(<*.-!)■ (20) 



a„ (21) 
f(oa-l), (22) 

/ , \ n-2 

n , , / I 



a 
a 



There will be a radiation-dominated primordial phase if q = 1/2. On the 
other hand, in the singular phase, the scale factor and its first derivative ap- 
proach constant values, and the second derivative, d, diverges as t — > t s , since 
n < 2. 

In terms of the conformal time, dn = a~ 1 dt, we have for the scale factor 
evolution to leading order: 

• Radiation-dominated primordial phase: 

a = aor). (24) 

• Singular phase: 

a = a s . (25) 
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The scale factor and its first derivative must be continuous during the tran- 
sition from one phase to the other. If a s is the scale factor value at the moment 
of the transition, and Hq the corresponding Hubble parameter value, then the 
transition moment is given by rj t — 1/{H a s ) and do = H a 2 s . 

The isotropic and homogeneous form we have assumed for the cosmological 
evolution of the scale factor, a(t), towards a sudden singularity captures the 
essential features of the general solution of the Einstein equations near such a 
singularity. For some more general scenarios, including anisotropies, see refs [211 

m m on 123 mi- 

We are ready now to determine the solution for the master equations for 
some relevant configurations. 



4.3 Big rip: massless, minimally coupled case 

When £ = m = it is more practical to solve directly the equation for the field 
<f>. The Klein-Gordon equation assumes the form, 

cb" + 2-0' + fc 2 = 0, (26) 

where from now on k denotes the wave number related to the Fourier decom- 
position of the spatial dependence of the scalar field. When the scale factor is 
given by a — ao|?7| 1+3a , representing the big rip scenario, the solution is: 

cj> k (r l ,x) = c l ri v H^\kr 1 y % - s , (27) 

where C\ is an integration constant, independent of k, v = — 2{i+zl) anc ^ Hi,(x) 
is the Hankel function of the first kind of order v. This form of the final solution, 
which should contain in principle two independent functions and coefficients 
with dependence on k, is dictated by the imposition of an initial Bunch-Davies 
vacuum state [27] . 



4.4 Sudden singularity: massless, minimally coupled case 

Let us return to the Klein-Gordon equation (|26[) . For the sudden singularity, we 
have previously considered two phases, a radiative phase and a singular phase. 
The initial conditions, corresponding to the Bunch-Davies vacuum state, are 
imposed at the beginning of the radiative phase, where the solution of the 
Klein-Gordon equation may be expresses in the form of plane waves: 

pikri 



But, the main interested is in the behavior of the quantum field near the singu- 
larity. In this phase, the solution of equation (f2l)]) reads, 

0(t,,2) = e -* \ A+e^-^ + A-e-^+k-^X ( 2 g) 
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where ui = \Jk 2 — 1 and the constants A± are fixed by the matching conditions 
across the two phases: 

A ± = l{l k e^ + \«±k). (30) 

4.5 Sudden singularity: massive, conformally coupled case 

If the mass is non-zero, and £ = i, the resulting Klein-Gordon equation takes 
the following form during the singular phase: 

X " + (fc 2 +m 2 a 2 ) X = 0. (31) 

where \ = a<p. During the radiative phase it is reasonable to consider the 
massless approximation, since we are in the regime where a — > 0. During the 
singular phase, the equation is 

x" + (k 2 +m 2 a 2 QV 2 )x = 0. (32) 

which again has plane wave solutions with a modified frequency. Hence, the 
solutions during the two phases can be written as follows: 

gifc?7 

4>k(v) = (primordial phase), (33) 

V 2k 



Mv) = Soie'*" + foe"* 5 " (singular phase), (34) 

where £01,02 are constants, to be fixed by the matching conditions, and uj = 
^/k 2 + m 2 ag. 



5 The energy of the created particles and the 
regularization procedure 

The — component of the energy-momentum tensor gives the energy density 
of a given configuration. The solutions shown in the previous section give the 
expression for the scalar field (massive or massless) for a given mode k. Inserting 
this expression in the energy-momentum tensor, the corresponding — com- 
ponent gives the energy associated with this mode. The total energy is obtained 
by integrating over all modes. In general, this integration leads to a divergent 
quantity, a common problem in quantum field theory. Hence, in order to make 
sense of the energy associated with the field, a renormalisation procedure must 
be employed. The renormalisation can be interpreted as a redefinition of the 
fundamental constants present in the original problem. In order for this pro- 
cedure to be physically meaningful, it is essential that the final result does not 
depend on how the renormalisation is performed. 

Hence, our problem here is to compute the expression 
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with pk = Tq. More precisely, we are interested in the final expression for the 
renormalized energy, p ren . It is this expression that will give the energy associ- 
ated with any particles created by the quantum processes near the singularity, 
and which may be relevant to the computation of any associated back-reaction 
phenomena in the evolution of the universe. 

In order to regularize the expression of the energy, we use the n-wave method 
exposed in the reference [28]. This method is based on the Pauli-Villars tech- 
nique used for quantum field theory in Minkowski space-time. First, let us write 
the energy as, 

y* oo 

p = I pk{k,m)k 2 dk. (36) 



Jo 

Let us define, 

p k = —pk(nk,nm), (37) 
n 

where n is a parameter that characterizes the order of the divergence. From 
this expression we construct the quantities, 

np 00 

I- 1 , t^-. (38) 



The expression for the regularised energy is given by, 

P r = Pk -El-El- l -El (39) 

where E® eliminates the logarithmic divergence, E k the quadratic divergence, 
and E 2 the quartic divergence - all those that are normally present in the 
energy-momentum tensor. This regularization of the energy corresponds to a 
full renormalisation of the coupling constants, as described in [29l [30] . 

The goal now is to determine the final expression for the regularized energy- 
momentum tensor, especially for the energy density that corresponds to the 
component — of this tensor. 

5.1 Big rip: the massless, minimally coupled case 

Using the solution (f2~T|) and the expression for the energy given by the — 
component of the energy momentum tensor (|15[) with m = 0, we obtain, 



Pk = Arf^x 2 



+ HW(x)HW(x)]. (40) 

The integration on all k modes reveals the existence of logarithmic, quadratic 
and quartic divergences. Using the n-wave regularization scheme, we write 



^ren 



poo 

/ x 2 (p k -E° k -El)dx 
Jo 



x 2 E l ° 9 dx. (41) 



k 
1/a 



11 



The last term is a modification of usual method in order to include the loga- 
rithmic divergence in the massless case, see [13]. The result is: 

where 1\ is a number. This expression describes how the energy of the created 
particles evolves with time. 

Taking the ratio between created particles and the background phantom 
fluid p x , we find 

n ren l ij-„ \ 

^octT 6 ^), (42) 

Px 

which diverges as the singularity at r\ = 0_ is approached. This may indi- 
cate that the evolution of the universe is modified by quantum effects and 
the singularity is avoided. But, in order to verify such effect, a full back- 
reaction analysis is necessary. Similar studies for the big rip singularity have 
been made in the context of quantum cosmologies, with inconclusive results, see 
references [23 GH 133 GH] ■ 



5.2 Sudden singularity for the massless, minimally cou- 
pled case 

Let us turn now to the sudden singularity case. As in the previous analysis for 
the big rip, we consider a massless, minimally coupled scalar field. In contrast 
to the big rip analysis above, an approximation considering the two phases is 
now necessary, and the initial vacuum condition is imposed in the first phase, 
which is that of a radiation-dominated universe. Using the solutions for the 
Klein-Gordon for the two phase given above, we find that the total energy is 
given by the following integral: 

p = Ae y J | (2k 2 - l)k 2 

— cosujy + wsinwj/j-, (43) 

where y — 2(1 — 77) and ui — y/k 2 — 1. The background constants are fixed such 
the singularity occurs at rj = 1, y = 0. 

Employing the n-wave regularization scheme, the regularized energy can be 
determined [15] : 

p™ = A e y{chi(-y) + ^\, (44) 

where Chi denotes the hyperbolic cosine integral function. The regularized 
energy decreases as the singularity is approached (y — > 0). Hence, the quantum 
effects are ineffective in preventing the singularity, at least for the massless scalar 
field case. 
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5.3 Sudden singularity for the massive, conformal cou- 
pling 

In this case, as for the previous one, two phases were considered. The initial 
vacuum state is fixed during the first phase, and the energy is computed during 
the second phase, in order to evaluate the possibility of a back-reaction effect 
on the evolution of the universe. 

Using the solutions (|33I34[) . the matching conditions and the expression for 
the energy-momentum tensor (| 14[) we obtain for the energy of the k mode: 



Pk 



1 



An integration over all /c-modes gives, 

/* oo 

p = / pkd 3 k = 7r 



fc 2 w( 1 - 



dk. 



(45) 



(46) 



This expression clearly diverges so it is necessary to regularise it. But, heuristi- 
cally, since it is a polynomial expression, it seems clear that after regularisation 
we must obtain zero. Hence, the particle production should not contribute to 
the energy-momentum tensor and the sudden singularity is unaffected by these 
quantum effects 

Note that the integral (|46|) admits an analytical solution: 



Pk d 3 k 



k 2 Cj \ 1 



dk 



w< k\J k 2 + m 5 



nr 
T 



k* 
2 



In 



2( k + \Jk 2 + m 



(47) 



with fh — moo. There is no infrared divergence, but there is a logarithmic 
divergence when k — > oo (ultraviolet limit). 
We have, 



Pk 



= \/k 2 + : 



It follows that 



(n) 
Pk 



2k 



pk- 



k 2 



Vk 2 



Hence, only the zero-order term survives, and leads to, 



Pk 



Pk 



El 



Pk- Pk = 0. 



(48) 
(49) 

(50) 



As we suspected, the renormalized energy is zero. There is no effect, and the 
quantum phenomena associated with the cosmological dynamics do not change 
the character of the sudden singularity or prevent its occurrence. 
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6 Conclusions 



It is known that quantum effects may play an important role near classical 
singularities, for example, for the big bang cosmological scenario. Even if the 
back reaction of these quantum effects on the classical evolution is still an open 
question, there are at least hints that quantum effects may lead to a dramatic 
deviation from the classical behavior [37] . 

Recently, new kinds of singularities have been identified in cosmology, the 
so-called future singularities. While the big bang singularity occurs in the origin 
of the universe, future singularities may mark the end of the universe. The 
interest in this kind of singularity has increased recently because of the unusual 
possibility that the presently observed accelerated expansion of the universe may 
be driven by a phantom fluid, which violates all energy conditions. However, it 
had been shown some time earlier that such future singularities may occur even 
if the energy conditions are not violated. 

The fate of quantum effects near future singularities has been reviewed here, 
with special attention to the big rip singularity and the sudden singularity. 
Concerning the big rip, the quantum effects may be relevant, but a more careful 
analysis of the back reaction process is necessary to decide under what conditions 
the big rip can be evaded by quantum effects [16] . In the case of the sudden 
singularity, the results obtained so far indicate that this singularity is robust 
against quantum effects because they are negligible in its vicinity. 

It must be remarked however that such studies has been carried out using 
special configurations of scalar fields. More general quantum fields must be 
analyzed and, of course, the problem of the back reaction must also be treated 
in greater generality. 
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